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AN  EQUATION  FOR  RAPID  CALCULATION  OF 
STAGNATION  POINT  RADIATIVE  HEAT  TRANSFER  (INCLUDING 
SHOCK  LAYER  RADIATIVE  COOLING  AND  NONGRAY  SELF-ABSORPTION) 

by 

John  D.  Anderson,  Jr. 

ABSTRACT:  A  closed-form  equation  is  derived  for  stagnation  point 
reentry  radiative  heat  transfer  accounting  for  the  combined  effects 
of  radiative  cooling  and  nongray  self-absorption  within  the  shock 
layer.  The  equation  can  be  applied  for  both  continuum  and  atomic 
line  radiation.  In  addition,  the  equation  is  shown  to  agree  favorably 
with  existing  numerical  data  for  stagnation  point,  continuum,  radia¬ 
tive  heat  transfer  for  a  wide  variety  of  conditions „  Also,  the 
equation  is  shown  to  apply  to  the  snd-wall  radiative  heat  transfer 
behind  a  strong  reflected  shock  wave  in  a  shock  tube.  Finally,  the 
equation  provides  a  rapid  means  of  obtaining,  by  hand,  reasonably 
accurate  engineering  estimates  for  reentry  radiative  heat  transfer 
including  shock  layer  radiative  cooling  and  nongray  self-absorption. 
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An  Equation  for  Rapid  Calculation  of  Stagnation  Point  Radiative  Heat 
Transfer  (Including  Shock  Layer  Radiative  Cooling  and  Nongray  Self- 
Absorption) 

This  report  presents  a  closed-form  equation  for  stagnation  point 
reentry  radiative  heat  transfer  accounting  for  the  combined  effects 
of  radiative  cooling  and  nongray  self-absorption  within  the  shock 
layer.  The  equation  can  he  applied  for  both  continuum  and  atomic 
line  radiation.  In  addition  the  equation  is  shown  to  agree  favorably 
with  existing  numerical  data  for  stagnation  point,  continuum,  radia¬ 
tive  heat  transfer  for  a  wide  variety  of  conditions.  Also,  the 
equation  is  shown  to  apply  to  the  end-wall  radiative  heat  transfer 
behind  a  strong  reflected  shock  wave  in  a  shock  tube. 

This  project  was  performed  for  Foundational  Research  under  FR-61. 


E.  F.  SCHREITER 
Captain,  jSN 
Commander 

L.  SCHINDEL 
By  direction 


NOLTR  68-56 


CONTENTS 

Page 

INTRODUCTION ....". . . . . .......  1 

ANALYSIS...... . 3 

Background . * . . . .  3 

Vacuum  Ultraviolet  Contribution. . 8 

Long  Wave  Length  Contribution... . 13 

Application  of  Formula . • . . . . .  16 

RESULTS . 18 

CONCLUSIONS .  23 

REFERENCES. . < . .  26 

TABLE  I 

Tabulated  Information  from  Engineering  Formula  and 
Comparison  with.  Numerical  Results .  21 


NOLTR  68-56 


ILLUSTRATIONS 

Figure  Title 

1  Effect  of  radiative  cooling  and  nongray  self-absorption  on 
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the  stagnation  region  shock  layer  enthalpy  profile. 

vw  -  50,000  ft/sec;  alt  -  200,000  feet 

4  Schematic  of  nongray  continuum  step  model  absorption 
coefficient. 
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Effect  of  radiative  cooling  and  nongray  self-absorption  on 
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and  B2  as  functions  of  temperature. 
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Comparison  of  engineering  formula  with  numerical  results 
for  stagnation  point  continuum  radiative  heat  transfer. 

Comparison  of  engineering  formula  with  measured  shock  tube 
end-wall  radiative  heat  transfer  rates  behind  a  reflected 
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NOMENCLATURE 


I , 


<0-1100 


Planck  black  body  function  -  2hc2£exp(hc/XkT)  -  ij  X"*® 

-1100A 

h 

J  ,B  dX 
J 1100 l  x 

total  radiative  energy  emitted  per  second  per  unit  volume 
evaluated  for  equilibrium  normal  shock  conditions 
radiative  energy  emitted  per  second  per  unit  volume  between 
0  and  1100A  and  evaluated  for  equilibrium  normal  shock 
conditions 

radiative  energy  emitted  per  second  per  unit  volume  above 
IIOOA  evaluated  for  equilibrium  normal  shock  conditions 
total  stagnation  point  radiative  heat  transfer 
stagnation  point  radiative  heat  transfer  below  1100A 


‘1100—* 


stagnation  point  radiative  heat  transfer  above  liooi 


nose  radius 

equilibrium  temperature  behind  a  normal  shock 
free  stream  velocity 

distance  measured  normal  to  the  surface  in  the  stagnation 
region 

shock  detachment  distance  for  the  shock  layer  without 
radiative  cooling 

exponential  integral  of  the  second  order;  Cn(t)  -  exp(-t/w)dw 

2  J0 

Pjf  s 
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F  radiation  loss  parameter,  Eq.  (6) 

Feff  effective  radiation  loss  parameter,  Eq,  (12) 

Tx  effective  radiation  loss  parameter  for  the  v-u-v  region 

Fg  effective  radiation  loss  parameter  for  the  long  wave  length 

region 

v-u-v  absorption  coefficient,  Eq.  (1) 
k2  long  wave  length  absorption  coefficient,  Eq.  (2) 

K.  spectral  absorption  coefficient 

X  wave  length 

p  density 

ry 

v-u-v  optical  length  -  J  K^(y)dy 

T^s  characteristic  v-u-v  optical  thickness  of  the  shock  layer  “ 

Kls  6A V 

Subscripts 

1  v-u-v  wave  length  region 

2  long  wave  length  region 

X  per  unit  wave  length 

«•  free  stream  conditions 

AD  without  radiative  coding 

s  equilibrium  conditions  immediately  behind  the  normal  bow 

shock  wave 
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INTRODUCTION 


Radiative  heat  transfer  from  the  high  temperature  shock  layer 
about  a  large,  blunt,  superorbital  reentry  vehicle  is  strongly  influ¬ 
enced  by  the  combined  effects  of  radiative  cooling  and  nongray  self- 

1-3 

absorption  within  the  shock-heated  gas.  In  fact,  recent  analyses 
have  shown  that,  for  superorbital  reentry  velocities,  radiative 
cooling  and  nongray  self-absoi’ption  can  reduce  the  stagnation  point 
radiative  heat  transfer  by  as  much  as  an  order  of  magnitude  in  com¬ 
parison  to  predictions  based  on  a  constant  property,  transparent  gas. 

Similar  effects  have  been  noted  on  the  high  temperature  gas  behind 

4  5 

a  strong  reflected  shock  wave.  *  Tf*»w  <wfrnir»«  l v  (n  nrHor  tn  t*k« 


Unfortunately,  in  order  to  take 


these  effects  into  account,  the  above  analyses  have  required  detailed 
and  te'Hou  numerical  calculations  of  the  radiating  shock  layer.  This 
situation  prompts  the  following  question:  can  a  simple,  approximate, 
closed-form  equation  be  derived  which  would  allow  rapid  but  accurate 
engineering  calculations  of  stagnation  point  radiative  heat  transfer 
taking  into  account  the  combined  effects  of  radiative  cooling  and 
nongray  self-absorption?  The  practical  benefit  of  such  an  engineering 
formula,  namely,  to  circumvent  lengthy  numerical  calculations,  is 
obvious.  In  answer  to  the  above  question,  the  oresent  paper  presents 
a  rational,  physical  derivation  of  such  a  closed-fora  equation  for 
stagnation  point  radiative  heat  transfer  accounting  for  radiative 
cooling  and  nongray  self -absorption.  This  equation  can  be  applied 
for  both  continuum  and  atomic  line  radiation.  In  addition,  the 
resulting  formula  is  shown  to  agree  favorably  with  existing  numerical 
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data  for  stagnation  point,  continuum,  radiative  heat  transfer  for  a 
wide  variety  of  conditions.  Also,  the  formula  is  shown  to  apply  to 
the  end-wall  radiative  heat  transfer  behind  a  strong  reflected  shock 
wave  in  a  shock  tube. 

The  present  analysis  is  an  outgrowth  of  previous  work  geared 
to  the  philosophy  of  simplifying  reentry  radiative  heat  transfer 
calculations  without  undue  sacrifice  in  accuracy.  In  particular, 
reference  (6)  documents  a  numerical  analysis  for  the  viscous,  nongray, 
radiating  stagnation  region  shock  layer,  using  approximate  step 
models  for  the  nongray,  continuum,  absorption  coefficient  of  high 
temperature  air.  In  fact,  the  explicit  purpose  of  this  analysis  was 
to  serve  as  a  numerical  instrument  with  which  to  investigate  the 
engineering  feasibility  of  such  step  model  absorption  coefficients. 
Early  results  ’  were  obtained  from  the  above  analysis,  using  a  now 
obsolete  step  model  absorption  coefficient.  More  recent  numerical 

g 

results  have  been  obtained  using  a  much  improved  two-step  model 
absorption  coefficient,  rationally  constructed  from  existing  quantum 
mechanical  data;  these  recent  results  have  shown  that  such  a  simple 
step  model  can  be  used  in  lieu  of  detailed  spectral  variations  in 
order  to  obtain  reasonable  engineering  results  for  shock  layer,  non- 

Q 

gray,  continuum  radiative  heat  transfer. 

The  closed-form  formula  for  reentry  radiative  heat  transfer 
derived  in  the  present  paper  represents  a  furtner  continuation  of 
this  engineering  philosophy.  The  formula  is  not  a  numerical  correla¬ 
tion  of  existing  data;  rather,  it  is  derived  on  a  rational  physical 
basis.  However,  a  hint  with  regard  to  a  crucial  physical  assumption 
is  revealed  by  close  examination  of  some  previously  unpublished  numeri¬ 
cal  results  obtained  during  the  preparation  of  reference  (8) ,  and  vvkes 
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a  simple  derivation  possible.  In  fact,  the  present  analysis  is  an 

illustration  of  one  role  computer  experiments  play  in  engineering 

9 

analyses,  as  described  in  a  recent  survey  by  Sichel. 

ANALYSIS 

Background 

As  background  for  the  following  derivation,  figures  1-3  show  , 

the  extent  to  which  radiative  cooling  and  nongray  self-absorption 
influence  the  radiating  stagnation  region  shock  layer  and  stagnation 
point  heat  transfer.  These  results  were  obtained  from  the  viscous, 
nongray,  continuum,  radiating  stagnation  region  analysis  which  is 

s 

mentioned  above,  and  which  is  simply  an  extension  of  an  earlier 

gray  gas  analysis  by  Howe  and  Viegas10  to  include  nongray  self-  s 

6  3  lo 

absorption.  These  analyses  are  well  documented,  ~  ’  u  and  therefore 

will  not  be  described  here.  It  is  sufficient  to  state  that  they 

contain:  (1)  coupling  of  the  radiative  energy  transport  with  the 

g&sdynamic  flow  field,  (2)  nongray  self-absorption  (in  the  case  of 

references  (6)  -  (3)),  (3)  a  fully  viscous  shock  layer  from  the  body 

to  the  bow  shock,  (4)  local  thermodynamic  and  chemical  equilibriur, 

and  (5)  a  self-similar  solution  limited  to  the  stagnation  region  of 

a  hypersonic,  thin,  radiating  shock  layer.  In  addition,  the  results 

in  figures  1-3  were  obtained  with  a  two-step  model,  nongray,  continuum 

absorption  coefficient  of  high  temperature  air  described  in  detail 

in  reference  (S) ;  these  results  have  been  shown  to  agree  favorably 

1  8 

with  the  detailed  spectral  calculations  of  Hoshizaki  and  Wilson.  * 

A  sketch  of  the  pertinent  step  model  absorption  coefficient  is  given 
in  figure  4,  where  and  K ^  are  the  vacuum  ultraviolet  (v-u-v)  and 

o 
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FIG.  1  EFFECT  OF  RADIATIVE  COOLING  AND  NONGRAY  SELF-ABSORPTION  ON 
STAGNATION  POINT  RADIATIVE  AND  CONVECTIVE  HEAT  TRANSFER. 

V  =  50,000  FT/SEC;  ALT  =  200,000  FEET 


.V#«fAv 
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* 


WITHOUT  RADIATIVE  COOLING 


y/  6 


flG.  3  EFFECT  OF  RADIATIVE  COOLING  AND  NONGRAY  SELF-ABSORPTION  ON 
THE  STAGNATION  REGION  SHOCK  LAYER  ENTHALPY  PROFILE. 

7^  =  50,000  FT/SEC;  ALT  =  200,000  FEET 
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is  obtained  from  the  results  of  Hahne,^*  and  Kg  is  obtained  from 

12 

one-half  the  radiance  values  of  Nar'”  't  al.  Reasonable  correla¬ 
tions  of  K^  and  k2  have  been  obtained  as  functions  of  the  local  gas 

g 

density  and  temperature. 

Figures  1-3  compare  results  obtained  with  and  without  radiative 
cooling  of  the  shock  layer,  i.e.,  with  and  without  coupling  of  the 
radiative  transport  with  the  gasdynamic  flow  field.  Nongray  self¬ 
absorption  is  included  in  both  cases.  These  results  illustrate  the 
strong  influence  of  radiative  cooling  and  nongray  self-absorption  on 

t 

radiative  heat  transfer  calculations;  consequently,  a  closed-form 
engineering  formula  for  super orbital,  reentry,  radiative  heat  transfer, 
Qr,  must  take  these  effects  into  account.  In  the  following  analysis, 
an  approximate,  closed-form  equation  for  QR  is  derived  which  indeed 
includes  the  combined  effects  of  radiative  cooling  and  nongray  self¬ 
absorption.  For  the  ease  of  continuum  radiation,  QR  will  be  considered 
as  the  sum  of  two  contributions:  (1)  the  v-u-v  radiative  heat  transfer 
below  HOOA,  Qro^xioo’  an<i  ^  th*  'l'ong  wave  length  radiative  heat 
transfer  above  1100 A,  QRxxqo-«' 
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Qr  “  QRq-1100  +  Qr1100-® 


(3) 


This  division  of  QR  into  two  distinct  (but  coupled)  parts  is  suggested 

by  the  spectral  variation  of  the  continuum  radiation  properties  for 

i  2 

high  temperature  air,  1  and  is  consistent  with  the  use  of  an  approxi- 

O 

mate  step  model  absorption  coefficient.  We  will  now  proceed  with 

the  derivation  of  closed-form  expressions  for  QRq_hqO  and  ^RllOO-“ 

which  allow  rapid  but  reasonably  accurate  calculations  of  these 

contributions,  Extension  to  the  case  of  atomic  line  radiation  will 

be  discussed  in  a  subsequent  section. 

Vacuum  Ultraviolet  Contribution 

Qr0-1100  is  raarkedly  affected  by  radiative  cooling  and  strong 

1—8 

self-abscrption  within  the  shock-heated  gas,  “  However,  these  two 
effects  combine  in  a  manner  that  allows  a  crucial  physical  assumption 
to  be  made,  subsequently  leading  to  a  simple  expression  for  Qr0_ixgo* 
To  illustrate  this  point,  figure  5  shows  the  effect  of  radiative 
cooling  and  nongray  self-absorption  on  the  local  values  of  the  product 
B.Go  as  *,  function  of  local  v-u-v  optical  length,  t,  ,  through  the 

1  *  f 1100 l  1 

shock  layer,  where  Bx  “  J  B^dX,  Gg  is  the  integro-exponential 

°  py 

function  of  second  order,  and  *  J  K^(y)dy.  The  areas  under  the 

curves  in  figure  5  are  proportional  to  the  v-u-v  radiative  heat 
transfer  to  the  surface  if  we  assume  one -dimensional  radiative  trans¬ 
fer  through,  the  shock  layer,  a  cold,  non-emitting,  black  surface,  and 

if  we  imply  that  the  gas  absorption  coefficient  is  adequately  repre- 

6  13 

sented  by  figure  4  (as  validated  in  reference  (8) ) .  That  is  * 

■*0-1100  -  2TJ*1S  VT>C2(T>dT  (4) 


8 


FIG.  4  SCHEMATIC  OF  NONGRAY  CONTINUUM  STEP  MODEL  ABSORPTION  COEFFICIENT 


(FT- LB) / (FT  -SEC-STER) 
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0-1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 


FIG.  5  EFFECT  OF  RADIATIVE  COOLING  AND  NONGRAY  SELF-ABSORPTION  ON  THE 
DISTRIBUTION  OF  €  2  THROUGH  THE  STAGNATION  REGION  SHOCK  LAYER 

AS  A  FUNCTION  OF  T, .  Vro  =  50,000  FT/SEC;  ALT  =  200,000  FEET;  R=2.0  FEET 
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The  previously  unpublished  numerical  results  shown  in  figure  5,  which 
were  obtained  during  the  preparation  of  reference  (8) ,  indicate  that 
the  local  maximum  of  is  predominantly  governed  by  radiative 

cooling,  whereas  the  decrease  in  for  larger  values  of  is  due 

to  self-absorption  (i.e.,  small  values  of  Cg) •  The  influence  of  the 
thermal  boundary  layer  near  the  wall  is  responsible  for  the  decrease 
in  BjCg  at  smaller  values  of  It  is  important  to  notice  in  figure 

5  that  radiative  cooling  and  self-absorption  flatten  the  distribution 
of  BjCg  through  the  shock  layer.  This  trend  is  further  supported  by 
the  results  of  figure  6,  which  graphically  shows  the  increasing  effect 


of  radiative  cooling  and  self-absorption  as  R  is  increased.  This 


flattening  of  the  B^Cg  curve  hints  very  strongly  at  a  simple  approxi¬ 
mation  for  Qr0_2iq0’  nainely>  tbe  area  under  the  curve  appears  to  be 
reasonably  approximated  by  a  rectangle  of  height  (B,C>)  and  of 

4.  4*  & 


length  Tlg, 


as  shown  by  the  dotted  lines  in  figure  6. 


(The  subscript 


s  implies  conditions  evaluated  immediately  behind  the  bow  shock.) 


With  this  approximation,  equation  (4)  becomes 


QRO-1100  “  2ir(BlC2)s  Tls  (5) 

This  expression  can  be  readily  evaluated  knowing  K^s  and  the  shock 
detachment  distance,,  <5,  as  will  be  outlined  in  a  subsequent  section, 
•In  addition,  the  results  shown  in  figure  6  indicate  that  the  above 
assumption,  and  therefore  equation  (5) ,  becomes  more  realistic  for 
increasing  values  of  R,  where  the  effects  of  radiative  cooling  and 
nongray  self-absorption  become  stronger,  and  where  the  thermal 
boundary  layer  becomes  a  smaller  fraction  of  the  total  shock  layer. 


(B^xlCf6  ( FT  -  LB )/  (FT2-  SEC  -  STEP ) 


‘vWVH  ' 


m 
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FiG.6  VARIATION  OF  € 2  WITH  THROUGH  THE  STAGNATION  REGION 
FOR  VARIOUS  NOSE  RADII .  V  =50,000  FT/SEC;  AIT=20G,000  FEET 
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Long  Wave  Length  Contribution 

In  contrast  to  the  strong  self-absorpt ion  occurring  in  the 
v-u~y  wave  length  range,  the  shock  layer  is  optically  thin  for  long 

•so 

wave  length  continuum  I'adi&tioo,  “  thus  providing  an  advantage  in 
obtaining  a  simple  expression  for  However,  Qr11g0  „  is 

affected  by  radiative  cooling  of  the  shock  layer  due  to  radiative 
emission  from  all  wave  length  ranges,  and  this  effect  must  be  taken 
into  account. 

For  a  transparent,  constant  property,  equilibrium  shock  layer,14 

the  stagnation  point  radiative  heat  transfer  is  given  by  E  6.^/2, 

s  AD' 

where  Eg  is  the  radiative  energy  emitted  by  the  shock-heated  gas 
per  unit  time  and  volume  evaluated  for  equilibrium  conditions  behind 
'the  normal  bow  shock,  and  6^  is  the  shock  detachment  distance.  How¬ 
ever,  if  the  transparent  shock  layer  is  cooled  due  to  radiative 
emission,  the  stagnation  point  radiative  heat  transfer  is  reduced,15 
and  is  given  by  (Es6AI/2)f (r) ,  where  f(D  si  and  is  a  function  of 
the  radiation  loss  parameter. 
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A  curve  fox*  f(D  has  been  obtained  by  Hoshizaki  (figure  9  of  referesc* 
(15));  as  indicated  in  reference  (6),  this  curve  is  apparently  indepej 
dent  of  the  gas  radlatiye  properties. 

From  the  above  experience,  it  appears  reasonable  that 
Qr1^00-“  oan  be  expressed  a h 
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where  E2g  *»  4s.'(x2B2)s  and  Bg  **  J  B^dX.  Thus, 


QRllOO-»  “  2ir(K2B2)s°AD  f(reff) 


*e£f  *s  “eriae{*  hare  as  an  "effective"  radiation  loss  parameter  which 
represents  the  ratio  of  the  radiative  energy  flux  out  of  a  constant 
property,  uncoupled,  nongray  shock  layer  to  the  flux  of  enthalpy 
ccnvected  into  the  shock  layer  from  the  free  stream  through  the  bow 
shock n  must  contain  the  influences  of  both  the  strongly  absorbed 

v-u-v  as  well  as  the  relatively  transparent  long  wavt.  length  ranges. 
Because  the  physical  effect  of  self-absorption  is  to  trap  some  of 
the  energy  within  the  absorbing  shock  layer  which  would  otherwise  be 
lost  if  the  gas  were  transparent,  reff  <  T.  Also,  because  of  the 
different  nature  of  the  two  wave  length  regions,  it  is  reasonable  to 
represent  reff  as  the  sum  of  two  terras 

reff  *  F1  +  r2  (9) 

where  and  r2  apply  to  the  v-u-v  and  long  wave  length  regions 
respectively.  An  expression  for  T2  is  immediately  suggested  by 
equation  (6)  applied  for  the  (transparent)  long  wave  length  region: 

r2  "  E2s6AD *  4*<K2B2>S6AD <10) 


An  expression  for  is  suggested  by  the  results  of  reference  (4) , 
which  indicate  that  an  "effective"  value  of  the  radiation  loss 
parasBeter  including  self-absorption  for  a  gray  gas  can  be  approxi¬ 
mated  by  the  product  £exp(-bTs> jr,  where  is  the  gray  optical 
thickness  of  the  shock  layer,  and  b  is  a  constant.  Applying  this 
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result  to  the  v-u-v  region  for  the  present  analysis,  and  approximating 
b  by  unity  as  suggested  by  inspection  of  the  numerical  results  of 
reference  (8),  the  following  expression  for  F^  is  obtained: 
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Note  that  F^  approaches  the  correct  optically  thin  expression  as 
Tls  'aecomQS  small.  However,  it  is  conceded  that  n  more  appropriate 
expression  for  F^  Might  be  obtained  with  further  study «  lu  fact, 
equation  (11)  shows  that  F,  approaches  zero  as  6._  (and  thus  t,  ) 
becomes  very  large,  whereas  it  is  physically  reasonable  chat  F-^ 
should  approach  a  finite  limit.  Consequently,  as  Tis  becomes  large, 
the  v-u-v  radiative  energy  less  through  the  bow  shock  is  not  entirely 
accounted  for  in  equation  (11) .  (On  tbs  other  hand,  it  will  b?  shown 
later  that  the  v-u-v  radiaxive  energy  loss  in  the  opposite  direction, 
towards  the  surface  c  does  indeed  approach  zero  as-  t^s  becomes  very 
large.)  Nevertheless,  for  practical  values  of  shock  layer  thickness, 
equation  (11)  is  considered  sufficient  for  the  present  purpose. 
Therefore,  from  equations  (9) ,  (10)  and  (11? , 
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Consequently,  equation  (8),  with  F^^  evaluated  as  equation  (12), 
allows  a  simple  evaluation  of  takes  into  account 

radiative  e&oling  fro«  all  wave  lengths. 
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with  r  given  by  equation  (12) .  For  a  given  reentry  trajectory 
point  (pw  and  VJ)  and  a  given  nose  radius -  R,  equation  (13)  can  be 
rapidly  evaluated  by  hand  in  a  straightforward  fashion  as  outlined 
below. 

(1)  Obtain  from  normal  shock  wave  tables  (such' as  reference 

(16))  the  equilibrium  density  and  temperature,  p  and  T  ,  behind  the 

s  s 

11 

bow  shock.  For  pg  and  T  ,  obtain  Klg  from  Hahnefs  tables  and  Kgg 

1? 

from  one  half  the  radiance  values  of  Nardone  et  al.  w  Alternatively, 
use  the  approximate  correlations  for  and  *23  available  from 
reference  (8) , 

(2)  Obtain  6^  from  hypersonic  flow  theory.  The  following 

17 

expression  is  recommended, 


6ad/r  - 


1  +  (2tV- 


where  e  -  Po/Ps* 

(3)  Calculate  Tlg;  vlg  -  ^is^AD*  Thsn  ^2^Tls^  from  existing 

iq  18 

tables  of  the  integro-exponential  function.  *'* 

(4)  For  T  ,  obtain  and  Bg  from  figure  7,  or  from  numerical 
evaluation  of  their  definitions. 


(5) 

(6) 


Evaluate  ?3if  from  equation  (12) . 


Obtain  f(reff)  from  figure  8, 


or  from  the  more  extended 


curve  in  figure  9  of  reference  (15) . 
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(7)  Evalute  the  stagnation  point  continuum,  nongray,  radiative 

heat  transfer  from  equation  (13) . 

The  above  development  has  been  carried  out  for  the  case  of 

continuum  radiation  with  the  continuum  absorption  coefficient  for 

high  temperature  air  reasonably  approximated  by  a  step  model  such 

as  that  shown  in  figure  4.  However,  this  analysis  can  be  extended 

to  include  atomic  line  radiation  as  reasoned  below.  Recent  measure-- 

ments*®”^  of  atomic  line  radiation  from  high  temperature  air  have 

indicated  that  the  line  intensities  are  strongly  self-absorbed.  In 

addition,  Nerem  is  currently  working  on  a  four-step  model  absorption 

coefficient  which  includes  atomic  line  as  well  as  continuum  contribu- 

22 

tions.  Also,  Olstad  at  NASA  Langley  is  developing  a  combined  atomic 
line  and  continuum  step  model  absorptiofa  coefficient.  Consequently, 
it  appears  that  equation  (13)  can  be  extended  in  a  straightforward 
manner  to  include  atomic  line  radiation  by  considering  a  multi-step 
absorption  coefficient  model  and  by  evaluating  the  contribution  to 
Qr  from  each  step  in  the  manner  prescribed  above.  Of  course,  Feff 
should  also  be  extended  to  account  for  each  step.  These  extensions 
are  not  explicitly  made  in  the  present  paper  due  to  the  current 
unavailability  of  a  precise  step  model  absorption  coefficient  which 
includes  atomic  line  radiation. 

RESULTS 

Equation  (13)  compares  favorably  with  existing  numerical  data 
for  stagnation  point,  nongray,  continuum,  radiative  heat  transfer, 
as  can  be  seen  from  figure  9  and  Table  1,  Figure  S  compares  the 
values  of  QR  obtained  from  equation  (13)  with  the  numerical  results 
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COMPARISON  OF  ENGINEERING  FORMULA  WITH  NUMERICAL  RESULTS 
FOR  STAGNATION  POINT  CONTINUUM  RADIATIVE  HEAT  TRANSFER 
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of  Hoshizaki  and  Wilson,  and  Anderson/  for  the  trajectory  point: 

Vw  *»  50,000  ft/sec;  altitude  =  200,000  feet.  Table  1  presents  a  more 
detailed  comparison  between  these  data,  including  additional  trajec¬ 
tory  points.  Even  though  the  existing  published  numerical  data  does 
not  extend  beyond  a  nose  radius  of  6  feet,  figure  9  shows  results  from 
equation  (13)  up  to  R  »  15  feet,  which  is  characteristic  of  Apollo 
size  reentry  vehicles.  These  comparisons  show  that  equation  (13)  is 
apparently  a  reasonable  engineering  formula  for  stagnation  point 
radiative  heat  transfer  taking  into  account  the  effects  of  radiative 
cooling  and  nongray  self -absorption  within  the  shock  layer.  In  addi¬ 
tion,  the  shape  of  the  curve  shown  in  figure  9  reflects  the  physical 
fact  that  Qr  is  dominated  by  the  v-u-v  contribution  at  small  R,  whereas 
for  large  R  the  v-u-v  radiation  is  strongly  self-absorbed  and  the  long 
wave  length  radiative  heat  transfer  is  by  far  the  dominant  contribu- 

1  4_« 

tion.  *  These  trends  can  be  seen  in  detail  in  Table  1.  In  fact, 
for  large  nose  radii  (R  on  the  order  of  15  feet)  the  present  approxi¬ 
mate  results  indicate  that  Qrq_iioo  may  be  on  the  order  of  one  Percent 
of  Q^.  These  proportions  remain  to  be  verified  by  detailed  numerical 
and/or  experimental  results  for  large  nose  radii;  however,  the  results 

5 

of  Nerem  and  Carlson  for  shock  tube  end-wall  radiative  heat  transfer 
behind  a  strong  reflected  shock  wave  in  air  indicate  that  such  propor¬ 
tions  between  Qr0_h00  and  ^R1100-»  are  reasonable  for  thick  shock 
layers.  Of  course,  on  a  physical  basis,  the  striking  reduction  of 
QRq  jqqq  for  large  nose  radii  is  due  to  the  strong  v-u-v  self-absorption, 
thus  attenuating  QRq_iioo*  and  to  radiativ8  cooling  of  the  shock 
layer,  which  reduces  the  local  radiative  emission  as  well  as  shifts 
the  peak  of  the  Planck  black  body  curve  to  longer  wave  lengths. 
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An  interesting  comparison  is  also  obtained  by  applying 

equation  (13)  to  the  end-wall  radiative  heat  transfer  behind  a 

strong  reflected  shock  wave  in  a  shock  tube;  for  this  case,  6^  - 

WRt,  where  is  the  reflected  shock  wave  velocity  and  t  is  the 

time  after  reflection.  Figure  10  compares  the  approximate  formula 

g 

with  the  shock  tube  results  of  Nerem  and  Golobic  for  the  case  of 
a  reflected  shock  wave  produced  by  an  initial  incident  shock 
velocity  of  8.85  mm/}*  sec  and  initial  driven  tube  pressure  of 
1  mm  Hg.  The  comparison  is  again  favorable,  even  though  in  this 
case  somewhat  dissimilar  quantities  are  being  compared.  That  is, 
equation  (13)  with  K^g  and  Kgg  evaluated  from  references  (11)  and 
(12)  respectively,  as  described  above,  represents  continuum,  nongray, 
radiative  heat  transfer  over  the  entire  wave  length  spectrum.  On 
the  other  hand,  the  shock  tube  results  of  Nerem  and  Golobic  give 
measured  radiative  heat  transfer  for  u.17  <  \  <  6p,  and  contain 
contributions  from  atomic  line  as  well  as  continuum  radiation. 

These  dissimilarities  are  to  some  degree  mutually  compensating, 
however,  the  relatively  favorable  agreement  shown  in  figure  10 
should  be  construed  as  somewhat  fortuitous.  Nevertheless,  the 
point  is  made  that  equation  (13)  can  be  used  to  predict  end-wall 
as  well  as  stagnation  point  radiative  heat  transfer. 

CONCLUSIONS 

An  approximate,  closed**form  equation  has  been  developed  which 
allows  rapid  calculation  of  reentry,  stagnation  point,  radiative 
heat  transfer  taking  into  account  the  effects  of  radiative  cooling 
and  nongray  self-absorption  within  the  shock  layer.  This  equation 
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FIG .  10  COMPARISON  OF  ENGINEERING  FORMULA  WITH  MEASURED  SHOCK 

TUBE  END-WALI.  RADIATIVE  HEAT  TRANSFER  RATES  BEHIND  A  REFLECTED 
SHOCK  WAVE 
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is  not  a  correlation  of  existing  numerical  data;  rather ,  it  has 
been  derived  on  a  rational  physical  basis.  However,  a  hint  with 
regard  to  a  crucial  physical  assumption  was  revealed  by  close 
examination  of  numerical  results;  namely,  that  the  product  3,€0 

A-  O 

is  reasonably  constant  through  the  snook  layer  due  to  the  combined 
influences  of  radiative  cooling  and  nongray  self-absorption.  The 
resulting  formula  compares  favorably  with  existing  numerical  calcu¬ 
lations  of  stagnation  point,  nongray,  continuum  radiative  heat 
transfer.  In  addition,  a  suggestion  is  made  for  an  extension  of 
the  formula  to  include  atomic  line  as  well  as  continuum  radiation, 
predicated  upon  the  future  development  of  a  step  model  absorption 
coefficient  which  represents  both  types  of  radiation.  Also,  in 
addition  to  the  stagnation  point  case,  the  formula  is  shown  to 
predict  end-wall  radiative  heat  transfer  behind  a  strong  reflected 
shock  wave.  Finally,  the  formula  provides  a  rapid  means  of  obtaining, 
by  hand,  reasonably  accurate  '•ngineering  estimates  of  reentry  radia- 
tive  heat  transfer  including  radiative  cooling  and  nongray  self¬ 
absorption,  thus  circumventing  lengthy  computer  solutions.  Of 
course,  the  formula  is  not  intended  to  replace  numerical  calculations 
when  very  high  accuracy  is  desired. 
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radiative  heat  transfer  for  a  wide  variety  of  conditions.  Also, 
the  equation  is  shown  to  apply  to  the  end-wall  radiative  heat 
transfer  behind  a  strong  reflected  shock  wave  in  a  shock  tube. 
Finally,  the  equation  provides  a  rapid  means  of  obtaining,  by  hand, 
reasonably  accurate  engineering  estimates  for  reentry  radiative  heat 
transfer  including  shock  layer  radiative  cooling  and  nongray  self¬ 
absorption. 


1.  radiative  heat  transfer 

2 .  stagnation  point 

3 .  nongray 

4.  reentry  heat  transfer 


INSTRUCTIONS 

of  «?eRS™G  A?TIVITY:  E<“«  name  and  addreaa  impose. 
finh  f  .f  it  ‘  ’  s“bcon‘rac‘°f.  grantee,  Department  of  De-  such  a: 
fenae  activity  or  other  organization  (corporate  author)  issuing 

report.  v*/ 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over- 
n<,,aV.fi,C“!l0n  the  ‘ePort-  Indicate  whether 
Restricted  Data’  is  included.  Marking  is  to  be  in  accord- 
ance  with  appropriate  security  regulations.  (3) 

26.  CJROUP:  Automatic  downgrading  is  specified  in  DoD  Di- 
rectlve  S200. 10  and  Armed  Forcer  Industrial  Manual.  Enter 
“^J,r°“P"umb*r-  Also  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author-  ^ 

*"*«, *he  comPiete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

U  a  meaningful  title  cannot  be  selected  without  classifica- 

tion,  show  title  classification  in  ail  capitals  in  parenthesis  (5) 

immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 

Give  tree  inclusive  dates  when  a  s^cuuc  reporting  period  is  c  V  * 
covered.  service 

5.  AITTHOR(S):  Enter  the  name(s)  of  authors)  as  shown  on  ,,  c„ 

°r  *be  teport*  Enter  last  name,  first  name,  middle  initial,  ,Y*  , 

If  military,  show  rank  and  branch  of  service.  The  name  of  tory  no 

the  principal  author  is  an  ahsolute  minimum  requirement.  12.  SP 

REPORT  DATE:  Enter  the  date  of  the  report  as  day,  *.*’'!  fe‘5, 

month,  year;  or  month,  year.  If  more  than  one  date  appears  ,n^  ‘or’ 

on  the  report,  use  date  of  publication.  13.  AB 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count  ?ummar 

should  follow  normal  pagination  procedures,  i.e.,  enter  the  *»  ”ay  ‘ 

number  of  priges  containing  information.  port.  If 

be  attar 

76.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report.  “  1 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter  an  indie 
the  applicable  number  of  the  contract  or  grant  under  which  formatic 

the  report  was  written. 

^  8d.  PROJECT  NUMBER:  Enter  the  appropriate  ever,  th 

military  department  identification,  such  aa  project  number, 
subproject  number,  system  numbers,  task  number,  etc.  H  KE 

9a.  ORIGINATOR’S  REPORT  NUMSER(S):  Enter  the  offi-  ^dc^cr 
cl •}  report  number  by  which  the  document  will  be  identified  selected 
snd  controlled  by  the  originating  activity.  This  number  must  fjers  su 
be  unique  to  this  report.  project 

96.  OTHER  REPORT  NUMBER(S):  If  the  report  has  been  words  b' 

assigned  any  other  report  numbers  (either  by  the  originetor  text.  T1 

or  by  the  sponsor )t  also  enter  this  numbers). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  “Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC.’’ 

(2)  “Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized. " 

(3)  **U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 

_  I* 

(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC  Other  qualified  users 
shall  request  through 

_ _ _  »» 

(5)  “All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 

_ _  II 

If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

1L  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay- 
trig  (or)  the  research  and  development.  Include  address, 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additional  space  13  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS).  (S).  (C),  or  (V). 

There  is  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

‘4  KEYWORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  modal  designation,  trade  ncme,  military 
project  code  name,  geographic  location,  may  be  used  as  key¬ 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rotes,  and  weights  is  optional. 


UNCLASSIFIED 

Security  Classification 


